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Abstract 

Using data collected by the fixed target Fermilab experiment FOCUS, we 
present several first measurements for the semileptonic decay D° — > K°ir~fi + is. 
Using a model that includes a K°n~ S-wave component, we measure the form 
factor ratios to be r v = 1.706 ± 0.677 ± 0.342 and r 2 = 0.912 ± 0.370 ± 0.104 
and the S-wave amplitude to be A = 0.347± 0.222 ± 0.053 GeV" 1 . Finally, we 
measure the vector semileptonic branching ratio T ^,^1 K = 0.337 ± 
0.034 ±0.013. 



1 Introduction 

Cabibbo allowed semileptonic decays have relatively large branching fractions 
and can be easily selected to achieve low levels of background contamination. 
The experimental results can be directly compared to theory where decay 
rates are calculated from first principles and include QCD effects in the form 
factors. Form factors are predicted in several models (quark models [1], Lat- 
tice QCD [2], and sum rules [3]). Once the form factors are determined, the 
CKM matrix elements can be calculated. While there have been several mea- 
surements of the D + [4,5,6,7,8,9] form factors, there are still no measurements 
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Fig. 1. Schematic of the decay D — > K*£ + u for angular variables definition. 



of the D° form factors for the vector semileptonic decays. We present the first 
measurement of the D° semileptonic form factor ratios for vector channels 
and the branching ratio T(D° -> K*(S92)-pL + u)/r(D° -> K°n-n + ). 1 Fur- 
thermore, we present an investigation of the S-wave component of the K°ti~ 
system and the measurement of its amplitude. This S-wave representation [10] 
was first used by FOCUS for the analysis of the decay D + — > K~n + fi + i/. 

The four-body decay amplitude can be parameterized by two masses and three 
angles. We use M(K°ir~), q 2 = (P M + P v ) 2 , and the three angles defined in 
Fig. 1: cosfly (the angle between the ir and the D in the K* rest frame), cos^ 
(the angle between the v and the D in the W rest frame), and x (the angle 
between the decay planes of the K* and the W). 

With these definitions the decay amplitude is written as: 2 

d 5 r 



dm KlT dq 2 d cos 9 v d cos 9 e dx 

(1 + cos0*) sm9 v e ix B K *-H + 
— (1 — cosOi) sm.Q v e~ lx B K *-H- 
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where K is the Kn system momentum in the D rest frame and B^*- and Ae zS 
are the Breit-Wigner and the S-wave components describing the spin one and 

1 Charge conjugation is implied throughout this letter. 

2 This model assumes that the q 2 dependence of the S-wave amplitude coupling 
to the virtual W + is the same as the Hq helicity amplitude describing the K* 
component. A study with as much as 100 times the statistics of this analysis has 
been performed for the FOCUS analysis of the decay D + — > K~tt + ji + u [10]. This 
study, where a significantly different kinematic dependence for the S-wave has been 
used, has reported a change in the form factors of less than 6% of the statistical 
error. 
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spin zero states of K°ir~, respectively. The four form factors (^-1,2,3 and V) 
that are necessary to describe a decay D — > V iv (where V stands for Vector), 
are included in the four helicity amplitudes: 

H±(q 2 ) = \{M D + m K MAq 2 ) =F J M °* V(q 2 )} (2) 

M D + TTlKn 

H (q 2 ) = 1 ^ [(M 2 - mL - <? 2 )(M D + m^A^q 2 ) 

AMI K 2 

° A 2 (q 2 )} (3) 



M D + m Kn 

tt , 2n M d K 2 M|, - rajL + g 2 , 

m K 7rV<l m d + m Kn 

+ M 2 / 4,(g 2 )] (4) 



A pole mass form is assumed for the form factors: 

= V(f) = (5, 



where M4 and My are the masses of the lowest cs states with the same quan- 
tum numbers as the W, namely M A = 2.5 GeV/c 2 and My = 2.1 GeV/c 2 
(which correspond to the masses of the and D* + , respectively). By in- 
cluding the parameter Al(0) in the constant that multiplies Eq. 1, the decay 
amplitude can be parameterized by the ratios of the form factors at q 2 = 0: 



v A 1 (oy 2 A^oy 3 Ax(o) 1 } 

We measure r v and r 2 . We have inadequate sensitivity to determine r 3 and 
we set its value to zero. From variations of this value we determine that the 
systematic uncertainty from setting r 3 = is negligible. 



2 Event Reconstruction and Selection 



FOCUS is a photoproduction experiment which collected data during the 
1996-1997 fixed-target run at Fermilab. The experiment, which is an upgrade 
of Fermilab experiment E687 [11,12], is characterized by excellent vertex res- 
olution and particle identification. For about 2/3 of the data taking the ex- 
perimental target was interleaved with a target silicon system [13] . The track 
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reconstruction downstream of the target is performed by four stations of sil- 
icon microstrips (SSD) and five stations of proportional wire chambers. The 
momentum of charged tracks is measured by the deflection in two magnets of 
opposite polarity. Charged particle identification is performed by three multi- 
cell threshold Cerenkov counters for electrons, pions, kaons, and protons [14]. 
Combining the information on the track momentum and the number of photo- 
electrons produced in the cells inside the (3=1 cone, a negative log-likelihood 
variable (W) for the hypothesis of the particle to be an electron, pion, kaon, or 
proton is determined. Particle identification is performed by a comparison of 
the probabilities for the different hypotheses and by requiring the hypothesis 
for the candidate particle to be higher than for the other hypotheses. Muons 
are identified by the hits left in tracking systems after penetrating approxi- 
mately 21 interaction lengths of shielding material [10]. 

In reconstructing D° — > K ir~fi + v, we select combinations of two charged 
tracks of opposite sign where one is identified as a pion and the other as a 
muon. For pion identification we require the pion hypothesis not to be dis- 
favored by more than six units of log-likelihood compared to the hypothesis 
with highest confidence level (min( W) — W w > —6), and to be favored by one 
unit of log-likelihood over the kaon hypothesis (to reduce the contamination 
from D° — > K~fi + v). For muon identification we require the track to have 
been reconstructed in the muon system (with at most one plane missed) with 
a confidence level greater than 1%. In order to reject background from the de- 
cays 7r + / ' K + — > fi + u, we require the muon trajectory to be consistent through 
the two analysis magnets with a confidence level greater than 1%. Each track 
must have momentum greater than 10 GeV/c. 

The two tracks are used to form the D° decay vertex, which is required to 
have C.L.> 5%, where C.L. is the confidence level. To reduce contamination 
from higher multiplicity decays, we require the probability for any other track 
reconstructed in the SSD system to come from the decay vertex to be lower 
than 0.1%. This requirement does not apply to the tracks used for the primary 
vertex reconstruction. To minimize background from hadronic re-interactions 
in the target, the decay vertex must lie at least one sigma outside of the 
target. The K° is reconstructed as a Kg from the decay K° s — > 7r~7r + [15]. The 
invariant mass is required to be within three sigma of the nominal Kg mass. If 
the pions are reconstructed using information from the silicon vertex detectors, 
the reconstructed Kg direction is used in the reconstruction of the D° vertex. 
In order to enhance the probability that our K s ir~ combination comes from 
a fT*(892)~, we require the reconstructed K s ir~ mass to be within one T of 
the nominal K*(S92)~ mass. The K*(S92)~ natural width T (50 MeV/c 2 ) 
is much larger than the experimental resolution on the reconstructed K s 7r~ 
mass (5 MeV/c 2 ). The invariant mass M(K s ir~ fi + ) is required to be lower 
than 1.8 GeV/c 2 . This cut significantly reduces combinatoric background since 
M(K s rc~ is kinematically limited to be below the nominal D° mass and 
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rejects D° — > Kg-n ir + decays when one of the pions is misidentified as a 
muon. 

We use the SSD tracks which have not been used in the D° decay reconstruc- 
tion to form primary vertex candidates. Each candidate is formed by starting 
with two tracks that make a vertex with C.L.> 1% and adding other tracks 
so long as the CL. remains greater than 1%. When a vertex is formed the re- 
maining tracks are used to form a second candidate in the same way and so on 
for the other candidates. We select the candidate with the highest multiplic- 
ity, and arbitrate ties by keeping the one with higher significance of separation 
from the secondary vertex. The significance of separation, which is given by 
the ratio of the distance between the two vertices divided by its error (£/ae), 
is required to satisfy ijan > 5. We "tag" the D° by requiring that it comes 
from the decay D* + — > D°ir+. The kinematics of this decay result in the pion 
having low momentum and being called a "slow" pion (tt s ). The pion must 
be one of the tracks used in the primary vertex reconstruction. It must have 
mm(W) - W w > -6 and p > 2 GeV/c. 



3 Fitting for the Form Factor Ratios and K°ty S-wave Amplitude 

In order to determine the form factor ratios r v and T2 we use a combined 
fit of the mass difference AM = M(D*) — M(D) and the three dimensional 
distribution cos 9 V vs. cos0£ vs. q 2 . For the AM component we use 60 bins 
in the region 0.14-0.20 GeV/c 2 . For the cos9 v vs. cos6>£ vs. q 2 distribution 
we select events with AM < 0.15 GeV/c 2 , and divide the phase space into 
four equally spaced bins for each of the two angular variables and two equally 
spaced bins for q 2 . The AM distribution, where signal and background events 
have a very different shape, is used to evaluate the background level. The 
binning choice for cosOy vs. cos0£ vs. q 2 gives information on the angular 
distributions of the W and the K°ir~ decays for two regions of q 2 . At low 
q 2 the angular dependence is more dramatic, while a more isotropic behavior 
is expected for high q 2 values, where the helicity amplitudes contribute with 
similar strength. 

Two methods are used to find the momentum of the missing neutrino. To com- 
pute the q 2 we use a ll D* cone" algorithm. By imposing energy and momentum 
conservation in the K*/i rest frame and by constraining the D and the D* to 
their nominal masses, the magnitude of p(D°) (which in this frame is equal to 
p{y)) is determined, but the direction lies on a cone. The direction is chosen 
by selecting the solution that gives the best x 2 when compared to the D° di- 
rection as given by the line connecting the two vertices. To compute the mass 
difference, we determine the neutrino momentum using the "neutrino closure" 
algorithm. This method is based on energy and momentum conservation for 
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the decay D° — > K tt~ and uses the nominal mass of the D° meson. The 
algorithm allows us to determine the neutrino momentum up to a two fold 
ambiguity, which is resolved by choosing the solution with lowest AM. Monte 
Carlo studies show that this choice is most often the correct solution. 

We use a binned maximum likelihood fitting technique with: 

£ = n 3 - , >< n c , (7) 

ijk b ijk- m °m- 

where (n^fe) is the number of observed (expected) events in the ijk th bin 
of the three dimensional distribution and S m (N m ) is the number of observed 
(expected) events in the AM distribution. The number of expected events is 
given by signal and background contributions. Non-charm backgrounds are 
essentially removed by the Ijoi requirement, by discarding events where the 
reconstructed decay vertex of the D° lies within one standard deviation from 
the target, and by the muon requirement. Contamination from charm decays 
is accounted for by using a Monte Carlo (that will be called MC B kg) which 
simulates all known charm decays other than our signal mode. The shapes for 
both distributions are taken from the distributions of the reconstructed events 
in MCbkg> an d their amplitudes are free to float. The background levels in 
the two distributions are tied by imposing that the yield of the MCbkg i n 
the cos 9 V vs. cos 9 f vs. q 2 distribution is equal to the area of the background 
shape in the AM distribution for AM < 0.15 GeV/c 2 . This corresponds to the 
selection cut imposed on the events in the cos#y vs. cos#£ vs. q 2 distribution. 
For the AM distribution, the signal shape is taken from Monte Carlo generated 
D° — > K°7i~fi + h / events. For cos 9 V vs. cos 9 e vs. q 2 the signal contribution to 
riijk is computed as the number of events generated in the bin ijk corrected 
by the efficiency for that bin. We calculate the generated number of event in 
bin ijk as a function of the fit parameters r v and r 2 using a Monte Carlo event 
weighting procedure based on Ref. [16]. For each Monte Carlo event generated 
in the bin ijk, we fill that bin with a weight given by the ratio of the decay 
amplitude in Eq. 1 evaluated for the fit parameters r v and r 2 over the decay 
amplitude evaluated for the input Monte Carlo values. 3 The signal yields in 
the cos 9y vs. cos 9p vs. q 2 and in the AM distributions are constrained in the 
same way as explained for the background events. 

The combined fit is shown in Fig. 2. The number of signal events is 175 ± 17. 
The x 2 P er degree of freedom in Fig. 2b is 32/27 which corresponds to a 



3 The FOCUS Monte Carlo simulation uses the D + — > K ir + fi + v form factor 
ratios and the S-wave parameters measured in [4]: r v = 1.504 ± 0.057 ± 0.039, 
r 2 = 0.875 ± 0.049 ± 0.064, A = 0.330 ± 0.022 ± 0.015, and 6 = 0.68 ± 0.07 ± 0.05. 
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Fig. 2. a) AM fit, b) cosOy vs. cos#£ vs. q 2 fit. Points with error bars are data, 
histogram is the fit, dashed line is the background component. 

confidence level of 22%. We measure the form factor ratios to be: 

r v = 1.706 ± 0.677 (8) 
r 2 = 0.912 ±0.370 (9) 



where the errors are statistical. 

The fit for the amplitude of the S-wave is performed with the same technique 
as the fit for the form factor ratios. We fix the form factor ratios to the values 
found above and, based on isospin symmetry, we fix the phase to 0.68, the 
value measured for the D + — > K~ir + /ji + v decay. As described in section 5, 
the possible bias due to this assumption is included in the systematic uncer- 
tainty evaluation. We find that A does not depend strongly on the phase. We 
measure: 

A = 0.347 ± 0.222 GeV -1 (10) 
where the error is statistical. 



4 The Branching Ratio T{D° -> K*(892)-(i + v)/T(D° -> K tt-tt + ) 



The branching ratio is measured by dividing the efficiency corrected yields 
of the two modes. The normalization mode D° — > K°7t~tt + is reconstructed 
following the same procedure and applying the same requirements as for the 



8 



1000 



800 



| 600 



400 



200 



(a 



YIELD=1918 ± 52 




1.14 0.16 0.18 0.2 

M(K jt~7tV)-M(K k~k + ) GeV/c 2 




M(K 7t"7tV)-M(K kk + ) GeV/c 



Fig. 3. AM fit for data (a) and Monte Carlo (b). The signal is fit to two Gaussian 
distributions, the background is fit to the threshold function in Eq. 11. 

D° — > K°7T~ fi + u mode (when possible), in order to minimize bias due to 
possible inaccuracies in the Monte Carlo evaluation of the efficiency for the Kg, 
which is reconstructed in a very different way from ordinary tracks. D* + tag, 
vertex reconstruction, Kg reconstruction, and particle identification (except 
muon identification) are the same as for the semileptonic mode. The n with 
the opposite charge of the 7i s must pass identical requirements as the n~ in 
the semileptonic mode. The trajectory of the n with the same charge as the 7r s 
must be consistent through the two analysis magnets, as we require for the 
In addition it is required to have min(W^) — > —6. For the hadronic mode 
we do not require the event to be in the mass window around the if* (892) ~ 
nominal mass. The invariant mass M(K s 7i^'K + ) must lie within 24 MeV/c 2 
of the fit D° mass, both for data and Monte Carlo. This window corresponds 
to a two sigma cut. The AM distribution is fit to two Gaussian distributions 
for the signal (in order to account for different resolutions) and the following 
threshold function for the background: 

BKG(AM) = a (AM - m^) 1 ' 2 + b (AM - m n ) 3 / 2 + c (AM - m n f' 2 (11) 



The fit for data and Monte Carlo are shown in Fig. 3. The yield from the fit 
to the data is 1918 ± 52 events. 

The efficiency corrected yield of D° — > K *(892)~fi + v is determined by correct- 
ing the efficiency corrected yield of D° — > K tt~ fi + u for the amount of K*~ 
in the K ir~ system. Since the Monte Carlo simulation uses the form factor 
ratios and S-wave parameters measured in the much higher statistics FOCUS 
analysis of D + — > K~n + ji~v, which are in excellent agreement with that mea- 
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sured in the present analysis, we estimate that the correction factors for the 
number of reconstructed events in data and Monte Carlo are the same, and 
therefore cancel out. The number of generated events in Monte Carlo (which 
is used to calculate the efficiency) must be corrected by the relative branching 
ratio T(D° -> K* (892)- n + v)/V{D° -> K°ir-/i + is) used for the Monte Carlo 
simulation. This can be evaluated by integrating over phase space the decay 
amplitude when the exclusive K* mode is generated, and dividing by the in- 
tegral over phase space of the decay amplitude for the inclusive mode. We 
calculate that in our Monte Carlo simulation the branching ratio is 0.95. We 
measure the branching ratio to be: 



where the error is statistical. 



5 Systematic Uncertainty Evaluation 

We carefully considered and evaluated many possible sources of systematic 
uncertainty in our results. Systematic bias can be generated by a poor Monte 
Carlo simulation of the detector performance, resulting in erroneous estima- 
tion of the efficiency. Also the particular choice for the fitting technique and 
parameters may bias the measurement. 

The accuracy and correct estimation of the errors reported by the fitting 
method is evaluated by repeating the measurement on a thousand samples 
obtained from fluctuating the bin entries in the data histogram. From the 
Gaussian distribution of the returned values for each measured quantity (form 
factor ratios, S-wave amplitude, signal and background yields), we conclude 
that the fit method is not affected by systematic bias and returns correct 
values for the errors. 

The Monte Carlo evaluation of the efficiency is investigated by repeating the 
measurements for different variations of the selection cuts. As expected, when 
the efficiency is correctly estimated (for our level of accuracy), the results are 
always stable within errors. We evaluate a possible bias due to the Monte 
Carlo simulation with the "split sample" technique, derived from the S-factor 
method used by the Particle Data Group [17]. The data is split into statisti- 
cally independent samples; for example, if the momentum simulation is being 
investigated, the data is split into distinct momentum regions. The measure- 
ment is performed on each sample for the observable x (e.g. r v ) and a x 2 for the 
hypothesis that the independent measurements are consistent is calculated. A 
poor consistency might result from a badly estimated efficiency with respect 




= 0.337 ±0.034 



(12) 
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to the momentum. We define poor consistency to be the case where % 2 > 1. 
In this case, the errors on the different measurements are scaled in order to 
return x 2 = 1, and we calculate a systematic uncertainty for the x measure- 
ment by subtracting in quadrature the statistical error from the scaled error 
on the weighted average of the independent measurements. Additional details 
are given in Ref. [18]. 

The bias from fitting choices is evaluated as the variance of measurements ob- 
tained by varying such choices. We vary the bin size both for the AM and the 
cos dy vs. cos Op vs. q 2 distributions. For AM, we also vary the fitting range. 
The r v and r 2 parameters are also evaluated setting the S-wave parameters 
to zero. The S-wave amplitude is evaluated for two additional values of the 
phase (at plus and minus one sigma from the reference value). For the r v , 
r 2 , and A fits, we include a variation on the fitting technique. This second 
fitting technique accounts for the efficiency in a different way. The efficiency is 
taken into account by using the weighting method on the reconstructed Monte 
Carlo events, instead of the generated events. For each event that passes all 
the selection cuts, the bin in which the event was generated in is filled with 
the weight described in Section 3. 

For the branching ratio measurement, we investigate the bias due to Monte 
Carlo input parameters by varying the form factor ratios and the S-wave 
values, and by varying the resonant structure of K°it~it + . Also, a less refined 
simulation of the hadronic trigger is investigated. The systematic bias from 
the model used in the Monte Carlo is evaluated as the variance of the three 
measurements found with these variations and the standard result. 

The total systematic uncertainty is given by the sum in quadrature of the 
uncertainties from the independent sources. Table 1 summarizes the results of 
the systematic uncertainty evaluation for all of the measurements. Including 
the systematic uncertainty we measure: 

r v = 1.706 ± 0.677 (stat) ± 0.342 (sys) (13) 
r 2 = 0.912 ± 0.370 (stat) ± 0.104 (sys) (14) 
A = 0.347 ± 0.222 (stat) ± 0.053 (sys) GeV" 1 (15) 



0.337 ± 0.034 (stat) ± 0.013 (sys) (16) 



6 Conclusions 



We have presented an analysis of the semileptonic decay D° — > K°ty 
using FOCUS data. Using a model which includes a K°n~ S-wave compo- 
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Table 1 

The systematic uncertainties from the Monte Carlo efficiency and acceptance eval- 
uation, the fitting condition, and total for r v , r 2 , A, and the branching ratio are 
shown. For the branching ratio, the systematic from the input parameters and trig- 
ger simulation in the Monte Carlo is also evaluated. 



Systematic Error 



Source 


c(r„) 


cr(r 2 ) 


ff(A)(GeV _1 ) 


cj(BR) 


MC Simulation 


0.0 


0.0 


0.0 


0.0 


Fit 


0.342 


0.104 


0.053 


0.002 


Model 








0.013 


Total 


0.342 


0.104 


0.053 


0.013 



nent that interferes with the dominant i^*(892) _ state, we have measured 
for the first time the D° form factor ratios for vector channels and the S- 
wave amplitude. We also report the first measurement of the branching ratio 
T(D ^K*(892)-h + is)/T(D ^K°it-it + ). 

Table 2 

The measurement of r v , r 2 , and A presented in this letter are compared to the 
FOCUS results for the decay D + — ► K~Tr + fi + u. We fix the S-wave phase to 0.68, 
the value measured for the D + . 







D+^K-tt+h+is 


T 2 

^(GeW 1 ) 


1.706 ± 0.677 ± 0.342 
0.912 ±0.370 ±0.104 
0.347 ± 0.222 ±0.053 


1.504 ±0.057 ±0.039 
0.875 ± 0.049 ± 0.064 
0.330 ± 0.022 ± 0.015 



From isospin symmetry, the expected values of r v , r 2 , and A can be directly 
compared to the results of the FOCUS measurements for the decay D + — > 
K~ir + fi + is, which uses the same model as the analysis presented in this letter. 
We find excellent agreement with the values for the D + , see Table 2. We 
calculate that in our model, where the K°ir~ system is given by a scalar and 
a vector component, the scalar fraction is 6%. 

The branching ratio value can also be estimated from the D + analysis using 
isospin symmetry: 

T{D°^K*-fi + u) _ t(D°) T{D+^ K*°fi + u) B(D+ ^K-tt+tt+) 
T(D ^K°7i-n+) ~t(D+) X T{D+^K-7i+7i+) * B(D° ^K°n-n+) ^ ' 

Since the decay dynamics do not depend on the lepton species, we compare the 
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Fig. 4. The T(D° -> K* (892)- ^+ 1/) /r (D° -> i^°7r-7r+) FOCUS measure- 
ment is compared to the CLEO-II measurement of the semielectronic mode 
r(D° -> K*(892)-e+i/)/r(D° -> K vr-7r + ), with the CLEO-c preliminary mea- 
surement of B(D° — ► X*(892)-e + ^ e ) divided by the Particle Data Group average 
for B(D° — > X°7r _ 7r + ), and to an estimate from isospin symmetry. The semielec- 
tronic results are corrected to account for the smaller electron mass when compared 
to the muon and they do not include the S-wave component. 

branching ratio result to measurements that use the semielectronic channel. 
Differences in the decay rate are only due to the larger mass of the muon 
as compared to the electron. In the semimuonic mode the phase space is 
reduced and there is a more significant contribution from the m 2 term of 
the decay amplitude (see Eq. 1). According to the PDG, the electron values 
should be corrected by a factor of 0.952 to compare to the muon results. We 
apply this correction and compare our results to the CLEO-II measurement 
of T(D° -> K*-e + v)/T(D° -> K°tt-tt + ) [19]. We also compare our results to 
the recent preliminary result from CLEO-c of the absolute branching fraction 
B(D° — > K*~e + u) (presented in conference proceedings [20]) divided by the 
PDG average of B(D° — > K°tt^tt + ). The comparison of our branching ratio 
measurement with the semielectronic results and with the calculation in Eq. 17 
is shown in Fig. 4. Only the calculation from isospin symmetry includes the 
effects of the S-wave component. The three estimates come from different 
measurements, and are in excellent agreement with each other and with our 
measurement. 
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